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Abstract. The temperature-concentration phase diagram of the Si-doped spin-Peierls compound CuGeOs 
is investigated by means of neutron scattering and muon spin rotation spectroscopy in order to determine 
the microscopic distribution of the magnetic and lattice dimerised regions as a function of doping. The 
analysis of the zero-field muon spectra has confirmed the spatial inhomogeneity of the staggered magneti- 
sation that characterises the antiferromagnetic superlattice peaks observed with neutrons. In addition, the 
variation of the macroscopic order parameter with doping can be understood by considering the evolution 
of the local magnetic moment as well as of the various regions contributing to the muon signal. 

PACS. 75.30.K Magnetic phase boundaries (including magnetic transitions, metamagnetism, etc.) - 
76.75.+i Muon spin rotation and relaxation - 75.25.+Z Spin arrangements in magnetically ordered mate- 
rials (including neutron and spin-polarised electron studies, synchrotron-source x-ray scattering, etc.) 



1 Introduction 

Since the first observation that the inorganic compound 
CuGeO-3 undergoes a spin-Peierls (SP) transition jjj, an 
extensive study of the doped system has been undertaken. 
Site (Cui-sM^GeOg) @,|,| and bond (CuGei^Si.Os) 
doping H[| studies have revealed the existence of a new 
antiferromagnetic phase below the SP ordering tempera- 
ture Tgp. Grenier et al. M have shown that the tempera- 
ture concentration(T-x) phase diagram for site and bond 
doping in CuGeOa exhibits a similar behaviour. The main 
features are the reduction of Tsp with increasing doping 
concentration, the suppression of the SP transition at a 
critical concentration x c , the onset of an AF phase char- 
acterised by an ordering temperature T^v increasing with 
x, and decreasing for concentrations larger than x c . The 
critical concentration x c is of the order of 3% for site dop- 
ing and about three times smaller for the bond doping 
scenario. 

An inspection of Mg-doped compounds near x c by 
means of magnetic susceptibility [0 and neutron scatter- 
ing measurements Q has provided a more detailed de- 
scription of the phase diagram for the site-doped SP. The 
ordering temperature Tn and the order parameter fJ- e ff 
of the antiferromagnetic phase exhibit a discontinuity at 
x c ?» 2.7%. The interpretation that has been proposed is 
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that doped CuGe03 undergoes a first order phase tran- 
sition at the critical concentration x c that constitutes a 
compositional phase boundary. A distinction is therefore 
introduced between a dimerised AF phase (D-AF) below 
x c , and a uniform AF phase for higher concentration (U- 
AF). Susceptib ility measurements near x c performed by 
Masuda et al. M] have revealed the presence of a dou- 
ble peak region, which was interpreted as a coexistence 
of both D-AF and U-AF phases. In addition they have 
performed high-resolution synchrotron diffraction on the 
(3/2 1 3/2) superlattice reflection. For low doping, the X- 
ray peak profile is resolution limited, indicating the pres- 
ence of long range order in the lattice dimerisation. In the 
region of doping which exhibits the double peak in the 
susceptibility, the superlattice X-ray peak broadens, thus 
indicating the onset of short range order in the SP phase. 

On the other hand, an extensive study of Si-doped 
compounds with susceptibility measurements was perfor- 
med by Grenier et al. JL1[ . They have shown that the 
introduction of Si ions reduces the intensity of the sig- 
nal attributed to the dimerised phase. Moreover an addi- 
tional peak appears below T$p, which is interpreted as 
the onset of an antiferromagnetic ordering at a tempera- 
ture Tjv- Above Tjy, they observe an additional intensity 
in the susceptibility which is attributed to the freeing of 
S=l/2 spins near the doping centers. The proportion of 
the SP signal linearly decreases with increasing doping, 



2 



F. Semadeni et al.: Macroscopic and Local Magnetic Moments in Si-doped CuGeOs 



whereas the proportion of free spins increases. However, 
no double peak feature is observed in the susceptibility 
measurements in the vicinity of the critical doping con- 
centration where the SP phase collapses. 

In order to explain the effect of impurities on the SP 
phase a theoretical model has been proposed by Saito and 
Fukuyama |l2],[lj|. A phase Hamiltonian that describes a 
onc-dimcnsional AF chain coupled with a 3D lattice dis- 
torsion field was used. The observation of both SP and 
AF signatures can be consequently explained by consid- 
ering a ground state where two long range order parame- 
ters coexist with a spatial variation. The lattice dimerisa- 
tion is minimal near the doping centers, where the stag- 
gered magnetisation takes its maximal value. Moreover, 
the phase transition from D-AF to U-AF is predicted to 
be of first order if the spin-phonon coupling is small when 
compared with the interchain interaction. A first experi- 
mental evidence of the spatial inhomogeneity for the mag- 
netic moments was seen with muon spin rotation experi- 
ments on Si- and Zn-doped CuGeC>3 compounds Q. 

In the present work we report a detailed study of the 
temperature-concentration (T-x) phase diagram of Si-doped 
CuGeC>3 single crystals by means of neutron diffraction 
as well as zero-field field muon spin rotation (/j,SR). The 
different signals observed with the neutrons, namely the 
AF and SP superlattice peaks, are interpreted in terms of 
volume fractions via the analysis of the zero-field muon 
spectra. Finally, the local magnetic moment measured by 
^tSR is compared with the macroscopic order parameter 
obtained by neutron diffraction. 



2 Experimental Details 

The Si-doped CuGeC>3 single crystals have been grown 
using the floating zone method. The impurity concentra- 
tion has been determined by means of inductively coupled 
plasma atomic emission spectroscopy (ICP-AES), with an 
accuracy of about 0.1%. The samples have been charac- 
terised by bulk susceptibility measurements which have 
been published elsewhere M . 

The elastic neutron scattering experiments have been 
performed on the three-axis spectrometer for cold neu- 
trons TASP, at the neutron spallation source SINQ, on 
a series of doped CuGei-^Si^Oa crystals (0.7% < x < 
3.8%). The incident neutron wave vector was kept fixed 
at ki = 2.662 A -1 . Higher-order neutrons have been sup- 
pressed by using a pyrolitic graphite filter. The samples 
have been oriented with the (Okl) zone in the scatter- 
ing plane, and mounted in an ILL-type cryostat, which 
achieves a base temperature of 1.5K. 

Zero field /iSR measurements have been performed on 
the spectrometers GPS and LTF, at the muon facility of 
PSI. The polarisation of the incident muon beam was par- 
allel to the a-axis of the sample. 



3 Results 
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Fig. 1. Profile for the magnetic superlattice peak (0 1 1/2) for 
the x=0.82%-doped sample, at 1.5K. A similar measurement 
at 5K shows the absence of higher order contamination. 



3.1 Neutron diffraction 

A typical profile of an antiferromagnetic superlattice peak 
is shown in Fig. [l] for the 0.82% sample. The tempera- 
ture dependence of the maximum intensity is presented in 
Fig. H In order to confirm the onset of the SP dimerisa- 
tion in the low doping regime, the reflection (1/2 3 1/2) 
was also measured. The magnetic intensity has been fitted 
with the following equation J3J: 
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which describes the critical behaviour of the magnetic in- 
tensity near the AF transition, weighted with a Gaussian 
distribution for T/v that accounts for inhomogeneities in 
the sample. N is a normalisation constant, (3 the critical 
exponent and BG the nonmagnetic background contri- 
bution. The exponent /3 remains almost constant over the 
whole diagram, with an average value of about 0.20±0.04, 
similar to that observed in Zn-doped crystals J3| . The inho- 
mogeneity in the concentration provides a AT^ of about 
0.3 ± 0.05K. The temperature dependence of the SP peak 
shown in Fig. hi was fitted with a function similar fel to 
the one of Eq. jj], 

The effective magnetic moments /x e / / have been calcu- 
lated by normalising the magnetic intensity at the satura- 
tion value with the nuclear structure factor of the (0 2 1) 
reflection, and corrected for the magnetic form factor of 
the free Cu 2+ ion |lq |. We emphasise that the neutron 
diffraction technique is a non local probe, and hence does 
provide a macroscopic order parameter, that is averaged 
over the whole volume of the sample. 

As a next step, selected samples (0.82%, 1.7%, 2.38% 
and 3.8%) belonging to the series of single crystals mea- 
sured with neutrons have been investigated by muon spin 
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Fig. 2. Temperature dependence of the (0 1 1/2) antiferro- 
magnetic and (1/2 3 1/2) spin-Peierls superlattice peaks in 
the 0.82% Si-doped CuGeOa single crystal. The solid line for 
(0 1 1/2) results from a fit according to Eq. A similar func- 
tion was used for the SP peak. 



rotation spectroscopy in the three temperature regions 
that have been determined with neutron and susceptibility 
measurements, according to the experimental T-x diagram 
presented in Fig. |^. The values obtained from the various 
techniques are found to be in good agreement, indicat- 
ing that the muon and the neutrons observe a magnetic 
ordering at the same transition temperature. The obser- 
vation of the SP dimerisation in susceptibility data for 
the samples in the low doping regime is also confirmed by 
neutrons. The discrepancy between both methods is due 
to the fact that the SP superlattice peak that is measured 
with neutrons becomes extremely small while approach- 
ing the critical concentration. Furthermore, we emphasise 
that Tsp is not observed above x—1.2%. 



3.2 ^SR results 

In the paramagnetic phase, the zero-field muon spectra 
have been fitted with the Kubo-Toyabe function which 
accounts for the muon depolarisation originating from the 
nuclear moments (the depolarisation rate is of the order 
of 0.08 MHz). By lowering the temperature, in the low 
doping regime, the muon spectrum can be still explained 
by the same relaxation function, as it was observed for 
pure CuGeC>3 |l7j]. In the magnetically ordered phase, be- 
low Tjv, the muon depolarisation P^t) in zero field was 
analysed over the whole doping regime with the following 
function (see also fllif): 



P»(t) = A r i x exp(-At) 

+ A osc cxp(-rt) • cos(27rz/i 
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Fig. 3. Temperature-concentration phase diagram of the Si- 
doped CuGeOs. The AF ordering temperature Tjv and the 
spin-Peierls transition temperature Tsp have been determined 
by bulk susceptibility jl5| and neutron diffraction experiments. 
The Tn determined by fj.SK measurements are also reported. 
The dashed lines are a guide to the eye. 



second term reflects a precessing part of the muon polari- 
sation, where r is the depolarisation rate and <f> the initial 
phase shift. The total observed amplitude (A r i x + A osc ) 
is significantly lower than the total amplitude observed 
above Tjv, indicating that part of the muon ensemble is 
depolarised within the dead time of the spectrometer (i.e. 
a depolarisation rate higher than ca. 150 MHz). 

According to the T-x diagram shown in Fig. ||, the 
emergence at low temperature of the precessing signal 
(A osc ) in zero-field is attributed to the ordering of the 
Cu 2+ moments near the doping centers. In this vein, the 
frequency is, which is proportional to the magnetic field 

at the muon site, mirrors the static ordering of the 
moments. 

The relaxing component in Eq. is ascribed to sam- 
ple regions where no coherent static magnetism is present. 
However, the similar temperature dependence of A and v 
(see Fig. ^ and the inset of Fig. p|) indicates that these 
regions have reduced dimensions, leading to a detectable 
influence of the ordered neighbouring regions on the field 
distribution Z\/7 M (where 7 M is the gyromagnetic ratio of 
the muon) at the muon sites. This supports the idea that 
there is no macroscopic phase separation in the samples 
over the whole doping range, similarly to what was ob- 
served by Kojima et al. for Si- and Zn-doped CuGeC>3 
compounds in the low doping regime |14 . Therefore, and 
similarly to the parameter v, the parameter A can be 
taken as a measure of the static ordered moments. Since 
A r i x 3> A osc the parameter A appears to be better deter- 



The first term in Eq. || describes a non-precessing part 
of the muon signal, which relaxes at a rate A, whereas the 



mined than v and consequently A will be utilised in the 
following discussion as a measure of the static moment. 

The doping dependence of the relaxation rate A, ex- 
trapolated at the saturation value, is shown in Fig. [| The 
dependence upon Si-doping is almost constant within the 
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Fig. 4. Temperature dependence of the relaxation rate A and 
the precessing frequency v for the 1.7% doped single crystal 
measured in zero-field /iSR (a). The magnetic moment mea- 
sured by means of neutron diffraction for the same sample is 
shown in (b). The solid lines are fits according to Eq. hi. 



precision of the experimental values. A exhibits however 
an increase above x=2.38%. According to the direct rela- 
tion between A and v, this behaviour indicates that within 
the doping range that has been investigated the local mag- 
netic field at the muon sites, and therefore the ordered 
moment increases only slightly with increasing tempera- 
ture. 

Figure || shows the evolution of the the normalised am- 
plitudes A r i x (light grey area) and A osc (white area), as a 
function of Si concentration. The amplitudes have been 
normalised with respect to the total muon asymmetry 
A to t, that was determined in the paramagnetic phase. We 
obtain that the amplitude corresponding to the magnet- 
ically ordered phase, A osc , represents about 10% of A tot . 
The relaxing part A r i x lies on the other hand between 40% 
and 60%. The part of A to t (dark grey area in Fig. ^) that 
could not be detected (« 40%) will be discussed below. 

In order to compare the local magnetic properties of 
the doped samples, as provided by the muons, with the 
macroscopic order parameter measured with neutrons, we 
use the quantity A ■ A osc . In this product, the local mag- 




0.00 0.01 0.02 0.03 0.04 0.05 
Si concentration x 

Fig. 5. Muon relaxation rate A extrapolated to T=0, plotted 
as a function of doping. The value is proportional to the local 
magnetic moment at the muon site (see inset). The dashed line 
is a guide to the eyes. The inset shows the ratio v/A for the 
x=1.7% sample, confirming the proportionality between both 
parameters (v/A « 2, as indicated by the dotted line). 
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Fig. 6. The normalised amplitudes of the muon depolarisation 
Arix (light grey area) and A oac (white area), obtained from 
Eq. as a function of Si doping. The quantities have been 
normalised to the total muon asymmetry A to t- The dark grey 
region refers to lost part of the muon signal (i.e. 1 - A r i x - 

A OS c) • 



netic moment (proportional to A) is hence weighted with 
the volume fraction (A osc ) of magnetised regions present 
in the sample. 

The dependence of the AF order parameter \i e $ t as a 
function of concentration is shown in Fig. [?], as obtained 
by both neutron diffraction and zero-field /iSR methods. 
Heff is seen to increase with the doping concentration, 
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Fig. 7. The order parameter fi e ff as a function of Si doping x 
measured with the neutron diffraction method (open circles). 
The effective magnetic moment determined by ^SR (explained 
in the text) is also shown (solid circles). The dashed line is a 
guide to the eyes. 



the spin chains, this distance corresponds to w 5 lattice 
units in both directions on each side of the muon site. The 
present measurements confirm the spatial inhomogeneity 
of the magnetic moments that was predicted theoretically 

60- 

The maximum in the order parameter [i e jf that was 
observed by neutron diffraction at x c =1.7% refers to a 
macroscopic value. In other words, the magnetic moment 
determined by this technique is averaged over the whole 
volume of the sample. The results provided by ^SR indi- 
cate that the volume fraction of the magnetically ordered 
regions as well as the boundary regions varies upon dop- 
ing. On the other hand, no anomaly is seen in the local 
magnetic moment near the critical value x c . We suggest 
therefore that the variation of the macroscopic AF order 
parameter upon doping is controlled by a variation in the 
volume fraction of magnetic domains rather than by a 
drastic change of the magnetic moments. 



5 Conclusion 



reaches a maximum at x c = 1.7%, where the SP signal is 
not observed anymore (see Fig. ||). 



4 Discussion 

The agreement between neutron and /zSR results (see Figs. || 
and [?]) allows to associate the precessing signal, referred 
to A osc , with the D-AF phase in the low doping regime, 
and with the U-AF phase in the high doping regime. 

The purely relaxing signal, A r i x , is attributed, for x < 
1.23%, to the regions where the lattice dimerisation is 
maximal. However, for higher doping concentrations, this 
amplitude does not disappear. It is thought to arise from 
an increase of static magnetic disorder that is induced by 
high doping. 

As the sum of A r i x and A osc constitutes not more than 
60% of the total muon asymmetry in the whole investi- 
gated doping regime, we may attribute the loss of signal 
to the boundary between the AF regions and the dimerised 
regions where a rapid spatial variation of the magnetisa- 
tion gives rise to large field gradient and therefore leads to 
a depolarisation rate too fast to be observed by the fiSR 
spectrometer. However, the neutrons are sensitive to these 
boundary regions that contribute to the intensity of the 
AF peaks. 

The /iSR measurements have shown on the other hand 
that the muon signal described by A r i x is influenced by the 
dipolar fields of the magnetically ordered regions and still 
undergoes a depolarisation. This is a clear sign that the 
domains contributing to the purely relaxing muon signal 
have a limited size. Considering the magnitude of the mag- 
netic moments that are determined by neutron diffraction, 
0.1 to 0.4 fiB, the maximal distance over which the muon 
still sees a dipolar field is estimated to be « 15 A. By con- 
sidering an homogeneous distribution of impurities along 



The compilation of information provided by neutron and 
/iSR have provided a detailed insight about the origin of 
the various magnetic and non magnetic contributions ob- 
served in the Si-doped CuGe03. The /iSR results have 
revealed that the regions that contribute to the antiferro- 
magnetic and lattice dimerisation peaks observed by neu- 
tron and susceptibility in the Si-doped CuGeOs have a 
typical size of 5 lattice units along the spin chains. 

Furthermore, the doping variation of the magnetic mo- 
ment measured with neutrons could be explained via the 
evolution of the volume fraction attributed to magneti- 
cally ordered regions near the impurity centers. 

One of the authors (F.S.) would like to thank the ETH- 
Council for the supply of a research grant. The U.S. De- 
partment of Energy is also acknowledged for financial sup- 
port (G.S.). 
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